THII  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOD  DIRECTIVE  5200,20  AND 
NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

r 

APPROVED  FOR  PUBLIC  RELEASE J 
DISTRIBUTION  UNLIMITED, 


Armed  Services  Technical  information  Agency 

Reproduced  by 

DOCUMENT  SERVICE  CENTER 

RNITT  BIILOIII,  MYTH,  2,  OHIO 


Thtt  document  la  lb*  property  of  the  United  States 
Government.  It  is  furnished  for  the  duration  o 4  the  contract  and 
shall  be  returned  when  no  longer  required,  or  upon 
recall  by  A8T1A  to  the  following  address: 

Ai  med  Services  Technical  Information  Agency,  Document  Service  Center, 

Knott  Building,  Dayton  2,  Ohio. 


NOTICE:  WHEW  GOVERNMENT  OR  OTHER  DRAWINGS,  SPECIFICATIONS  OR  OTHER  DATA 
AkU  USED  FOR  ANT  PURPOSE  OTHER  THAN  IN  CONNECTION  WTTH  A  DEFINITELY  RELATED 
GOVERNMENT  PROCUREMENT  OPERATION,  THE  U.  8.  GOVERNMENT  THEREBY  INCURS 
NO  RESPONSIBILITY,  NOR  ANY  OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE 
GOVERNMENT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE 
SAID  DRAWD  G6,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  AMY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 
PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE, 
USE  OR  SELL  ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


MU 


HI! .  ''MM*' 


FC 


il  II II II II II II II  ii 


REPORT  No.  962 


Some  Aspects  Of 
Three-Dimensional  Boundary 
Layer  Flows 


R.  SEDNEY 


DEPARTMENT  OF  THE  ARMY  PROJECT  No.  5B03-03-001 
ORDNANCE  RESEARCH  AND  DEVELOPMENT  PROJECT  No.  TB3-O108& 

BALLISTIC  RESEARCH  LABORATORIES 


ABERDEEN  PROVING  GROUND.  MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


REPORT  NO.  9&2 
September  1955 


SOME  ASPECTS  OP  THREE-DIMENSIONAL  BOUNDARY 
LAYER  FLOWS 

R.  Sedney 


Department  of  the  Army  Project  No.  5B03-03-001 
Ordnance  Research  and  Development  Project  No.  TB3-0108 


ABERDEEN  PROVING  GROUND, 


MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


REPORT  NO.  962 


RSedney/mjf 

Aberdeen  Proving  Ground,  Mi. 
September  1955 


SOME  ASPECTS  OF  THREE-DIMENSIONAL  BOUNDARY  LAYER  FLOWS 

ABSTRACT 

The  equations  for  laminar  boundary  layer  flow  over  a  general  surface 
in  three-dimensions  are  analyzed  in  a  normal  coordinate  system.  The 
invariance  properties  of  these  equations  are  found  using  the  concept  of 
subtensors.  The  boundary  layer  equations  are  not  tensor  equations  but 
subtensor  equations.  Conditions  for  the  Cartesian  form  of  the  equations 
sure  given  and  a  criterion  for  no  secondary  flow  is  found  in  terms  of  the 
geodesics  of  the  body  surface.  The  displacement  effect  of  the  boundary 
layer  la  also  discussed. 
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SYMBOIS 


a. 


ap 


G1'G2 


H. 


k 

n 

or 


12  3 

u  ,  u  ,  u^. 


\>  Ug,  My 


or  u,  v,  v 
U,  V 


or 


x,  y,  z 


6  1,6  2 

x  *  x 


Surface  metric  tensor 
Space  metric  tensor 

Components  of  vector  whose  surface  divergence 
is  zero 

Surface  metric  for  orthogonal  coordinates 

Space  metric  tensor  for  orthogonal  coordinates 

Geodesic  curvature  of  surface  streamlines 
Unit  normal  vector  to  a  curve 

Boundary  layer  velocity  components 


Tangential  components  of  external  flow  velocity 

Normal  coordinates 
"Mass  flow  defects” 


A 

X 

v 

X 


Displacement  surface 

Unit  tangent  vector  to  a  curve 

Kinematic  viscosity 

Ratio  of  pressure  to  (constant)  density 


Subscripts  and  superscripts 

i,J,  etc.  Have  the  range  1,  2,  3. 

a,p,  etc.  Have  the  range  1,  2. 
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INTRODUCTION 


In  recent  times  interest  in  three-dimensional  boundary- layer  flows, 
measured  by  the  number  of  published  papers  on  the  subject,  has  grown 
considerably.  A  review  of  the  subject  is  given  by  Sears[l]  .  Most  of 
the  work  cited  by  Sears  is  concerned  with  the  solutions  to  particular 
problems.  To  the  author’s  knowledge  the  only  general  discussions  of 
boundary  layers  in  three-dimensional  flows  are  contained  in  the  work  of 
Howarth  [2]  ,  Moore  [3]  ,  and  Hayes  [u]  .  Hovever  there  Is  also  some  work 
of  0.  C.  Lin,  contained  In  Chapter  18  of  j>]  ,  m  which  the  appropriate 
equations  are  derived  without  much  discussion.  This  latter  work  seems 
to  be  the  first  to  derive  the  boundary  layer  equations  for  a  general  curved 
surface  but  has  been  overlooked  by  most  authors.  The  way  in  which  Lin 
derives  the  equations  turns  out  to  be  useful  in  discussing  the  invariance 
properties  of  the  boundary  layer  equations. 

In  this  paper  it  is  first  shown  that  Lin’s  approach  carries  over  to 
a  more  general  class  of  coordinate  systems.  Certain  invariance  properties 
follow  from  examining  the  resulting  equations.  Now  it  is  known  that  the 
two-dimensional  boundary  layer  equations  are  not  tensor  equations  and  one 
would  expect  the  same  result  in  three -dimensions  (in  [4j  ,  Hayes  seems  to 
imply  the  contrary).  However,  in  the  language  of  [6]  ,  the  boundary  layer 
equations  are  subtensor  equations  which  means  that  they  are  invariant  with 
respect  to  certain  types  of  coordinate  transformations.  This  does  not 
contradict  the  work  of  Lagerstrom  and  Kaplun,  [7]and[8j  ,  where  the  non¬ 
tensor  character  of  the  boundary  layer  equations  is  used  to  define  an 
"optimum  coordinate  system" .  A  similar  procedure  could  be  developed  for 
three-dimensional  flows,  however  a  more  general  coordinate  system  than 
that  used  here  would  have  to  be  considered. 

From  the  subtensor  form  of  the  equations,  it  is  easily  shown  that 
they  reduce  to  the  "Cartesian  form"  (as  is  true  for  any  two-dimensional 
flow)  for  any  surface  whose  Gaussian  curvature  is  zero  if  appropriate 
coordinates  are  used.  Howarth,  [2]  ,  concludes  that  this  happens  only  for 
planes  and  cylinders.  The  curvature  effects  that  Howarth  describes  include 
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property  that  the  distance  expression  Is  given  by 


'2  =  +  g33 ^to3)2 


ds  = 


that  is, 


*3a*° 


There  is  no  unique  way  of  choosing  our  normal  coordinate  system  since  all 
that  is  required  is  that  the  given  surface  he  one  of  the  family.  A  simple 
special  case  is  the  geodesic  normal  coordinate  system  in  which  the  family 
of  surfaces  is  obtained  by  measuring  off  constant  distances  from  a  given 
surface  along  the  geodesics  which  cut  the  surface  orthogonally.  In  this 
special  case 


This  is  the  coordinate  system  used  by  T.-tn  , 

The  surface  x^  ■  0  is,  in  general,  a  Riemannian  2- space.  In  this 
subepace  we  can  carry  out  tensor  operations  and  these  will  be  intimately 
related  to  the  tensor  operations  of  the  parent  3- space.  We  use  a  comma 
to  denote  covariant  differentiation  in  the  3- space  and  a  semi-colon  for 
the  same  operation  in  the  subspace.  A  set  of  quantities  Ta  which  trans¬ 
form  according  to  .the  law 

Tto  =  T0  (d*e/dx'a) 


under  a  transformation  of  the  form 


*•“  .  f(a>  (x*5) 


»  x' 


is  called  a  subtensor.  If  a  tensor  T^  is  split  up  into  two  groups  T^,  T^ 
then  for  (l)  TQ  is  a  subtensor  and  T^  is  a  subinvariant,  and  similarly 
for  higher  order  tensors.  Also  the  Christoffel  symbols,  in  which  one  or 
more  of  the  indices  have  the  value  3,  are  subtensors.  A  more  detailed 
discussion  of  subtensors  is  given  in  [6j  . 
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the  effects  of  a  curvilinear  cooruiaate  system,  (in  [3]  ,  Moore  states 

that  the  Cartesian  form  applies  to  any  curved  surface  which  is  certainly 
not  true . )  . 

Using  the  streamline  coordinates  as  in  [  4]  ,  it  follows  that  a 
boundary  layer  flow  is  essentially  two-dimensional  if  there  is  no 
secondary  flow.  A  simple  criterion  for  no  secondary  flow  is  given  in 
terms  of  the  geodesics  of  the  surface.  Finally  the  displacement  effect  of 
the  boundary  layer  is  discussed.  This  is  done  in  a  different  way  than 

by  Moore,  £9]  >  and  some  differences  between  the  two-  and  three-  dimensional, 
cases  are  discussed.  \ 

Only  laminar  flows  of  an  incompressible,  non-conducting  fluid  are 
discussed.  Some  of  the  conclusions  should  hold  in  more  general  cases, 
however.  Also  body  force*  have  been  neglected  and  difficulties  due  to 
separation  effects  are  not  discussed. 

THE  BOUNDARY  LAYER  EQUATIONS 

In  applying  the  boundary  layer  concept  to  the  flow  of  a  viscous 
fluid  we  have  same  given  surface  to  consider,  for  example  a  solid  body 
or  an  interface  between  tvo  fluids.  This  basic  surface  is  used  to  define 
a  coordinate  system:  it  is  one  of  the  coordinate  surfaces.  A  convenient 
class  of  coordinate  systems  is  the  normal  coordinate  systems,  [6]  ,  which 
is  defined  as  follows.  A  given  one-parameter  family  of  surfaces  in  a 
Riemannian  N- space  has  a  family  of  orthogonal  trajectories  such  that, 
undtfr  very  general  conditions,  only  one  trajectory  passes  through  each 
point  of  space.  The  parameter  of  the  family  of  surfaces  is  denoted  by  xN 
and  on  one  of  the  surfaces  a  coordinate  system  x1,  ...,  xN_1  is  set  up. 

We  shall  be  concerned  with  a  Euclidean  3- space  although  the  derivation  of 
the  boundary  layer  equations  would  proceed  in  the  same  way  for  more 
general  spaces. 

The  convention  is  adopted  that  Greek  suffixes  have  the  range  1,  2  and 
small  Latin  suffixes  the  range  1,  2,  3.  The  given  surface  over  which  we 
consider  the  boundary  layer  flow  is  the  one  on  which  the  coordinates  xa  are 
defined  and  for  this  surface  x^  =  0.  Normal  coordinate  systems  have  the 
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expression  is  given  by 
+  g33( dx3)2 


There  is  no  unique  way  of  choosing  our  normal  coordinate  system  since  all 
that  is  required  is  that  the  given  surface  he  one  of  the  family.  A  simple 
special  case  is  the  geodesic  normal  coordinate  system  in  which  the  family 
of  surfaces  is  obtained  by  measuring  off  constant  distances  from  a  given 
surface  along  the  geodesics  which  cut  the  surface  orthogonally.  In  this 
special  case 

®33  *  1 

This  is  the  coordinate  system  used  by  Lin  . 

The  surface  x3  =  0  is,  in  general,  a  Riemannian  2- space.  In  this 
subspace  we  can  carry  out  tensor  operations  and  these  will  be  intimately 
related  to  the  tensor  operations  of  the  parent  3- space.  We  use  a  comma 
to  denote  covariant  differentiation  in  the  3- space  and  a  semi-colon  for 
the  same  operation  in  the  subspace.  A  set  of  quantities  TQ  which  trans¬ 
form  according  to  .the  law 

l'a  =  Tp  (d^/^a) 


property  that  the  distance 

,2  ,  a.  6 

ds  =  dx  dxp 


that  is. 


=  0 


under  a  transformation  of  the  form 
x'a  -  f<°>  <**>) 
x'5  = 


(1) 


is  called  a  subtensor.  If  a  tensor  T^  is  split  up  into  two  groups  T^,  T^ 
then  for  (l)  T^  is  a  subtensor  and  is  a  subinvariant,  and  similarly 
for  higher  order  tensors.  Also  the  Christoffel  symbols,  in  which  one  or 
more  of  the  indices  have  the  value  3,  are  subtensors.  A  more  detailed 
discussion  of  subtensors  is  given  in  . 
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The  derivation  of  the  boundary  layer  equations  is  straight-forward 
but  tedious.  The  procedure  here  is  essentially  the  same  as  Lin's  for  the 
special  case  of  a  geodesic  normal  coordinate  system  end  therefore  will 
not  be  discussed  in  detail,  (dome  errors  were  found  in  the  details  of 
£5^  but  these  only  affected  terms  which  drop  out  in  the  boundary  layer 
approximation. )  Starting  from  the  Navier-Stokes  equations  in  tensor  form 


(duj/dt)  +  uJui;j  *  vgJkUi,Jk 


(2) 


where  s  is  the  ratio  of  pressure  to  the  constant  density,  the  following 
steps  are  necessary.  The  momentum  equations  are  split  up  into  two  groups 
(as  illustrated  above  with  T^);  the  covariant  derivatives  are  expressed 
in  terms  of  the  "sub- covariant  derivatives"!  a  transformation 


x3/vl/2 


V 


.1/2 


(3) 


is  applied;  all  quantities,  including  the  metric  tensor  and  Christoff  el 
symbols,  are  expanded  in  a  power  series  in  v1^2;  the  equations  for  the 
lowest  order  terms  then  yield  the  boundary  layer  equations 

(dUg/dt)  +  uPua.p  +  (Uj/gjjXdji^C)  *  -  *.a  +  (&2V^2)/g33 

*  0  (4) 

uP.p  +  CpuP  +  (dUj/dO/gjj  -  0 


where 


ca  "  (de55/dxa)/2g55 


and  the  metric  tensor  for  (4)  is  where 

^  fit 

00 

The  boundary  layer  equations  (4)  reduce  to  those  given  by  Lin  for  g^j  *  1. 
(Note  that,  in  the  procedure  outlined  above,  nothing  is  implied  about 
the  higher  approximations  obtained  from  the  series  expansion.  Consideration 
of  these  involves  significant  difficulties,  [lj]  .) 
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Under  transformations  of  the  form  (l)  it  is  easily  seen  that  equations 
(4)  are  invariant  since  and  g^  are  invariant .  Thus  the  boundary  layer 
equations  have  subtensor  character.  This  is  not  too  surprising  since  vhat 
destroys  the  tensor  character  of  this  approximation  is  the  transformation 
(3)  and  this  is  not  affected  by  (l).  Also  in  the  subtensor  form  it  becomes 
obvious  that  the  boundary  layer  equations  reduce  to  "Cartesian  form"  for 
any  surface  of  zero  curvature  (flat  2- space)  in  which  the  coordinates  x 
are  Cartesian  coordinates  provided  that  a  geodesic  normal  coordinate  system 
is  used,  i.e.  »  1.  This  conclusion  differs  from  Knwarth's,  [V]  , 
because  he  restricts  his  analysis  (which  uses  a  geodesic  normal  system)  to 
coordinate  systems  xa  which  can  be  Cartesian  only  for  planes  and  cylinders. 

Because  of  the  subtensor  character  of  (4)  there  is  complete  freedom  in  the 

rv  a 

choice  of  x  .  However,  the  choice  of  Cartesian  coordinates  x  on  a 

developable  surface  may  not  always  be  the  most  advantageous.  For  example, 

the  form  of  the  terms  x .  which  are  determined  by  the  external  flow,  may 

be  complicated  by  such  a  choice .  It  must  be  remembered  that  the  conclusion 

concerning  when  the  boundary  layer  equations  reduce  to  Cartesian  form 

applies  only  to  within  the  approximations  of  the  standard  theory.  For 

example,  for  the  flow  over  on  open-ended  cylinder  the  equations  are  in 

Cartesian  form(for  appropriate  xa)  but  the  usual  boundary  layer  approximations 

become  invalid  as  the  distance  from  the  leading  edge  increases.  Finally,  it 

may  be  noted,  even  for  flow  over  a  plane  the  equations  will  have  curvature 

terms  appearing  if  xa  are  not  Cartesian  coordinates. 

The  boundary  layer  equations  (4)  can  be  written  in  the  more  conventional 
form  in  terms  of  the  physical  components  of  the  velocity.  Before  doing  this 
let  us  specialize  to  orthogonal  coordinates  x  and  change  the  notation  for 
the  metric  tensor 


=  hQ2  (no  summation)  (5) 

gjjU^x^O)  -  h52 
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vhere  h  is  not  a  function  of  x^.  Also  denote  the  physical  components 

a 

of  the  velocity  by  u,  v,  and  w.  Thus 
u  *  ui/hi 

v  -  Ug/hg 

w  -  Uj/hj 

since  for  (4)  the  metric  tensor  components  are  the  lower  case  h's.  In 
terms  of  the  physical  components,  using  x,  y,  and  z  as  coordinates,  equations 
(4)  become 

2 

Ut  +  uu^/t^  +  vu^/hg  +  uvlt^y/l^hg  -  V  h^/h^hg 

+  wuz/h3  -  -  +  vuZ8  A i 

vt  +  *  vv/h2  +  uvhac/hlh2  -  u2‘Vhlh2 

+  w^hj  -  -  *  wM  A?  (6) 

*«  *° 

+  IWV]  ^  +  V*/hJ  ’  0 

vhere  the  subscripts  t,  x,  y,  and  z  denote  partial  differentiation!  These 
reduce  to  Hovarth's  equations  [2J  for  hj  ■  1,  i.e.  a  geodesic  normal  coordinate 
system.  Hsyes,  [V]  ,  gives  these  equations  (if  compressibility  is  neglected 
in  his  equations)  except  that  he,  in  some  wav,  allows  hj  to  depend  on  s. 

Finally,  it  can  be  remarked  that,  to  examine  the  flow  in  the  neighbor* 
hood  of  a  stagnation  point,  as  Howarth  has  done,  |lo]  ,  it  is  only  necessary 
to  Introduce  a  Riemannlan  coordinate  system  for  x°  with  origin  at  the 
stagnation  point, 

SECONDARY  FLOW 

Since  the  boundary  layer  equations  are  invariant  under  transformations 
of  the  surface  coordinates  x°,  it  is  natural  to  look  for  coordinates  which, 
under  certain  conditions,  simplify  the  general  equations.  An  ingenious 
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choice  of  coordinates  was  made  by  Hayes,  [VJ  .  For  a  steady  external 

flow  the  coordinate  is  chosen  along  the  streamlines  of  the  external 

o 

flow  evaluated  on  the  surface  and  x  is  along  the  orthogonal  trajectories 
of  these  streamlines.  If  U  and  V  are  the  physical  components  of  the 
external  flow,  evaluated  on  the  surface,  in  the  x  and  y  (or  x1  and  x2) 
directions,  then  for  this  choice  of  coordinates  V  ■  0  and 


-  «A  ■  WA 

-  *  A  -  - A 


(7) 


where 

k  -  Va 

i.e.  k  is  the  geodesic  curvature  of  the  streamlines.  Consider  steady  flow 
in  the  boundary  layer.  The  boundary  conditions  for  v  are 
v  =  0,  x  — *»oo 

end,  if  the  surface  over  which  the  flow  is  considered  is  a  solid,  non- 
spinning  body, 

v  =  0,  z  ■  0 

Now  if  k  *  0,  it  is  seen  from  (7)  and  the  second  equation  of  (6)  that 
v  -  0  is  a  solution.  The  term  "secondary  flow  '(also  cross  flow)  is  used 
to  indicate  that  the  streamlines  in  the  boundary  layer  flow  do  not  coincide 
with  the  external  flow  streamlines  evaluated  on  the  body.  Also  k  «  0  has  a 
simple  geometric  interpretation:  the  streamlines  are  geodesics  of 
the  surface.  Thus,  for  steady  flow  over  a  non- spinning  body,  there  is  no 
secondary  flow  if  the  external  flow  streamlines  are  geodesics  of  the  body 
surface. 


A  simple  example  of  the  above  conclusion  is  provided  by  the  flow  over 
a  flat  plate  with  an  arbitrary  leading  edge  placed  in  a  uniform  stream  at 
no  angle  of  attack;  see  Figure  1. 


FIGURE  I 
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The  external  flow  streamlines  are  Just  straight  lines  on  the  plate  so 
there  is  no  secondary  flow.  (Moore,  £$]  ,  reaches  this  conclusion  by 
finding  a  solution  of  the  flow  equations.  He  also  indicates  how 
difficulties  arise  if  the  leading  edge  does  not  have  a  continuously 
turning  tangent.)  A  special  case  of  this  is  the  yawed  flat  plate. 

However  consider  yawed  Infinite  cylinders,  which  many  authors  have  done, 
DJ.  It  is  easy  to  see,  using  the  above  criterion,  that  the  yawed  flat 
plate  is  the  only  case  with  no  secondary  flow  since  for  any  other  cylinder 
the  external  flow  streamlines  cannot  be  geodesics. 

Since  v  *  0  in  the  streamline  coordinates  for  k  ■  O,  the  boundary 
layer  flow  is  essentially  two-dimensional.  Equations  (6),  written  for 
a  geodesic  normal,  system,  become 

uVh!  +  +  vuzz 


(h2u)x/hih2  +  vz  •  0 


A  new  coordinate,  X,  can  be  introduced, 
x 

X  ■  f  h1  dx 

Jo 


and  then  equations  (8)  are  seen  to  be  in  exactly  the  form  of  the  boundary 
layer  equations  over  a  body  of  revolution.  Therefore  the  same 
transformation  that  Mangier,  [llj  ,  has  introduced  will  reduce  equations 
(8)  to  the  standard  two-dimensional  equations.  Mangier4 s  transformation 


(9) 


u  ■  u’ 


w 


V  -  b2J\ 
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Note  that  the  coordinate  transformation  here  is  not  of  the  type  (l) 
and  that  the  transformation  to  new  velocity  components  (u',w* )  does 
not  follow  a  tensor  law.  Hayes  has  Indicated  the  possibility  of 
transforming  the  no-secondary  flow  equations  to  two-dimensional 
form.  However  his  transformation  is  made  by  choosing  the  form  of  the 
metric  tensor  in  a  suitable  manner  and  it  would  be  very  difficult 
to  find  out  exactly  what  the  transformed  coordinates  are.  From  the 
discussion  above  it  is  seen  that  the  transformation  is  just  that  of 
Mangier,  except  that  y  appears  as  a  parameter. 

DISPLACEMENT  SURFACE 


The  displacement  of  the  external  flow  streamlines  by  the  retarding 

action  of  the  boundary  layer  is  an  important  effect.  In  two-dimensional 

flow  the  definition  of  displacement  thickness,  which  is  a  measure  of 

this  effect,  is  straightforward  and  there  are  several  equivalent 

definitions.  In  three-dimensional  flow,  if  one  proceeds  in  strict 

analogy  to  the  two-dimensional  case,  two  displacement  thicknesses 

and  &  0  where 
x*' 


can  be  defined,  f?]  .  These  are  8^ 


U8  «  fh  (u  -  u)  dx5 
ar  •  o 

(10) 

V5x2  *X  (V  -  v)  dx5 


16 


in  which  (u,v)  and  (U,V)  are  the  physical  components  of  the  velocity  in 
12 

the.  x  and  x  directions  in  the  boundary  layer  and  the  external  flow 
respectively  and  h  is  "some  location  veil  outside  the  boundary  layer." 

(The  coordinates  xa  are  again  general  coordinates,  not  the  streamline 
coordinates  of  the  previous  section. )  Moore  refers  to  these  lengths  as 
characterizing  mass  flow  defects  and  then,  using  these,  defines  a 
displacement  surface  as  follows.  In  a  geodesic  normal  coordinate  system, 
the  displacement  surface  x^  *  A(x^,x2)  is  an  "impermeable  surface  which 
would  deflect  a  nonviscous  fluid  in  such  a  way  as  to  produce  a  normal 
velocity  (W)  satisfying" 

W  ■  w(x^,x2)  at  x5  *  h(x*,x2) 

where  h  has  the  same  meaning  as  above,  and  W  is  the  external  flow  normal 
velocity.  After  some  approximations  Moore  gives  a  differential  equation 
for  A,  [9]  . 

Here,  using  a  slight  modification  of  one  of  the  two-dimensional 
definitions,  an  equation  for  A  is  obtained.  This  approach  is  quite 
different  from  Moore's  but,  making  Just  the  boundary  layer  approximations, 
this  equation  reduces  to  that  of  Moore.  Analysis  of  the  original  equation 
shows  an  interesting  difference  between  the  three-  end  two-dimensional  cases. 

In  the  following  derivation  geodesic  normal  coordinates  cure  used  and 

all  vectors  are  in  physical  components.  It  is  convenient  to  use  the  suffix 

3 

notation  for  the  range  1,  2  but  we  set  z  ■  x  .  First  we  dispose  of  some 
geometrical  preliminaries.  It  can  be  shown  that  (see  ,  for  example) 
the  metric  tensor  components  g^  are  quadratic  functions  of  z 

-  “op  +  2W  +  (11) 

where  a^  has  already  been  defined  and 

*<#  =  <8V'8*)o  /2 

°af>  *  <A*i/aiV2' 


IT 


Also  the  following  relation  holds 


c  -  m  a0Db  bQ 

ap  ao 


(12) 


Specializing  to  orthogonal  surface  coordinates ,  by  means  of  (12), 
equations  (ll)  can  be  written  as  a  perfect  square.  In  the  notation 
of  (5) 


H. 


a 


h  +  i  z 
a  a 


(13) 


where 


*  bg^/h^  (no  summation) 


(14) 


If  on  a  surface  z  =  c  =  constant  we  draw  a  simply  chosed  curve,  the 
tangent  vector,  X,  and  the  normal  vector,  n,  have  the  cosqponents 

X:  (HjdxVds,  Hgdx^ds) 

(15) 

n:  (-Hgdx^/ds,  HjdxVds) 


where  S  Is  the  arc  length  along  the  curve. 

To  define  a  displacement  surface  we  consider  the  flux  of  fluid 
through  a  developable  surface  S  formed  by  the  normals  to  the  surface  z  »  0, 
passing  through  a  simply  closed  curved  KQ  on  z  »  0>  see  Pig.  2. 
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FIGURE  2 
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This  flux  is  computed  for  that  part  of  S  between  x  =  0  and  z  ■  h.  For  a 
velocity  distribution  Q  this  flux  Is  given  by 


r® "  { 4'nds  ■  Cf  <wjtl ■ 


)dz 


(16) 


where  K  Is  the  trace  of  the  developable  surface  on  the  surface  z 
c 

o<c<h.  The  displacement  surface 


c  and 


z  ■  A(x\x2) 


Is  a  surface  such  that  the  flux  through  S,  for  0<z<h,  Is  the  same  for 
the  following  two  velocity  distributions 


Is 


II: 


*L 

*L 

<2 

*2 


P  ■  0  for  0<z<A 
U 

for  A<z<h 
V 

u 

for  0<z<h 
v 


From  (16)  we  set 


and  this  is  the  condition  imposed  to  find  A.  Interchanging  the  order  of 
integration  and  making  use  of  (13)  and  (14),  after  which  the  line  integrals 


are  calculated  for  Kq,  this  requirement  yields 


K 


(G^dx2  -  Ggtylx1)  -  0 


(17) 


where  ^ 

Gx  =  UA  [  1  +  (i2 A/2hg)  ]  +  f  (u  -  U)  [l  +  (zi g /d^)]  dz 

°h  ,  (lQ) 

Gg  =  VA  [l  +  (^A/21^)]  +  J  (v  -  V)  [l  +  (zi^)]  dz: 

The  integral  in  (17)  can  be  written  as  an  integral  over  the  area  bounded  by 
Kq  using  the  "surface  divergence  theorem"  [12]  .  Since  KQ  is  arbitrary 
we  then  obtain  the  result 


div  G  =  0 
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(19) 


for  the  vector  0: 
divergence'.'  £l2j  , 

div  0 


(0l,02)  where  the  operator  div  is  the  "surface 

1  §6  e 


(“iV]  V2 


or  in  the  subtensor  notation  (I9)can  be  written 


For  given  velocity  distributions  I  and  II  and  a  value  of  h,  (l<?)  is  a 
differential  equation  for  A.  However ,  in  boundary  layer  theory,  h  cannot 
be  given  a  definite  value.  In  fact,  in  accordance  with  the  requirements 
of  this  theory,  we  must  let  h-*©**’  (Extrapolating  from  the  two-dimensional 
case,  the  integrals  should  converge  since  u-?U,  v->V  exponentially  as 


Thus  far  we  have  made  no  approximations.  We  take  (U,V)  to  be  the 

1  p 

(physical)  velocity  components,  in  the  x  and  x  directions,  of  the  external 
flow  evaluated  on  the  surface  z  »  0  and  (u,v)  the  corresponding  velocity 
components  throughout  the  boundary  layer.  Applying  a  transformation  of 
the  type  (3)  to  z  and  A  and  keeping  only  the  lowest  order  terms  in  (l8) 
gives  the  boundary  layer  approximation  to  and  Gg 

G1  =  UA  +  C*  (u  -  U)dz  «  U(A  -  6  ^ 

J  _  /aaN 


Gp  =  VA  +  f  (v  -  V)dz  =  V(A  -  B  2) 

Jo  * 

using  the  definitions  (10).  With  these  expressions  (20)  for  the  vector  G, 
(19)  is  essentially  Moored  equation, fojt  for  the  displacement  surface. 
Moore  gives  some  examples  of  the  computation  of  A  for  special  kinds  of 
flow.  Here  we  give  a  simple  example  by  specializing  only  the  coordinate 
system.  For  the  streamline  coordinates  of  the  previous  section  V  *  o, 
therefore  Gg  vdz  and  (19)  becomes 


therefore 


d  hgU(A  -  /by}  =  -  d^Gg  )/dx2  =  F(x1,x2) 
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In  us 

h2U(A  -  8  )  =  /F(x1,x2)dx1  +  f(x2) 

2  * 

vhere  f(x  )  is  a  constant  of  the  integration  with  respect  to  x1  and  must 
be  evaluated  from  the  conditions  of  the  flow.  Note  that  in  the  streamline 

coordinates  8  2  is  undefined, 
x 

The  method  used  above  to  obtain  A  is  a  modification  of  the 
corresponding  two-dimensional  definition  ofi  displacement  thickness.  In 
two  dimensions  the  developable  surface  S  is  taken  to  be  a  plane 
perpendicular  to  the  plane  of  flow  and  as  a  result  no  differential 
equation  need  be  solved  for  the  displacement  thickness.  The  A 
determined  from  (19)  using  (20)  differs  -  from  the  6  of  (10)  by  a 
constant  of  integration  which,  in  most  cases  of  interest,  is  aero. 

This  is  discussed  by  Moore,  . 

However,  it  is  interesting  to  compare  the  general  expressions  (l8) 
for  the  vector  G  in  the  two-  and  three -dimensional  cases.  Let  x2  ■  0 
be  the  plane  of  flow  for  a  two-dimensional  flow.  Then  v  -  V  «  0  and 
it  is  easy  to  show  that  ig  »  0.  Thus  the  terms  that  are  neglected  after 
making  the  boundary  layer  approximations  to  obtain  (20)  disappear  auto¬ 
matically  for  the  special  case  of  two-dimensional  flow.  This  would  be 
important  if  one  wanted  to  calculate  approximations  to  a  viscous  flow 
beyond  the  classical  boundary  layer  theory  as  Kuo,  [l}J  ,  has  done  for 
two-dimensional  flow  past  a  flat  plate.  For  two-dimensional  flow  the 
displacement  thickness  expression  does  not  change  for  the  higher  order 
approximations  whereas  for  three  dimensional  flow,  including  axially 
symmetric  flow,  an  expansion  of  the  expressions  (18)  in  powers  of  v1/2 
would  lie  necessary. 

R.  SEDNEY  ' 
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